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ABSTRACT 

The formation of spectroscopic binaries (SB) may be a natural byproduct of star formation. The 

_4 . early dynamical evolution of multiple stellar systems after the initial fragmentation of molecular 

CL|. 

clouds leaves characteristic imprints on the properties of young, multiple stars. The discovery 



and the characterization of the youngest SB will allow us to infer the mechanisms and timescales 



involved in their formation. Our work aims to find spectroscopic companions around young stellar 

objects (YSO). We present a near-IR high-resolution (R ~ 60000) multi-epoch radial velocity 

survey of 7 YSO in the star forming region (SFR) p Ophiuchus. The radial velocities of each 

source were derived using a two-dimensional cross-correlation function, using the zero-point 

established by the Earth's atmosphere as reference. More than 14 spectral lines in the CO Av 
""NT ■ 

£f) . = (0-2) bandhead window were used in the cross-correlation against LTE atmospheric models 

to compute the final results. We found that the spectra of the protostars in our sample agree 
well with the predicted stellar photospheric profiles, indicating that the radial velocities derived 
are indeed of stellar nature. Three of the targets analyzed exhibit large radial velocity variations 
during the three observation epochs. These objects - pending further confirmation and orbital 
characteristics - may become the first evidence for proto-spectroscopic binaries, and will provide 
important constraints on their formation. Our preliminary binary fraction (BF) of ~71% (when 
merging our results with those of previous studies) is in line with the notion that multiplicity is 
very high at young ages and therefore a byproduct of star formation. 

Key words, infrared: stars - stars: formation - stars: pre-main sequence - binaries: close - 
spectroscopy 



1. Introduction 

Since the early twenties that there is growing evidence that multiple stars ma y be the rule and not 



the exception. This conception is more recentlt supported by observations of lDuquennov & Mavoi 



* Based on observations collected with the CRIRES spectrograph at the VLT/UT1 8.2-m Antu Telescope 
(ESO run ID. 081.C-0395(A)) at the Paranal Observatory, Chile 
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(1199 ll) who found that the stellar multiplicity among main sequence (MS) stars in the solar neigh- 
borhood may be as high as ~ 60%. Since then, more refined surveys probed the dependencies of the 
properties of stellar systems on physical parameters like age, mass, separation, and environment. 
These observations are highly relevant, because there is still no comprehens ive theory that explains 
the properties of multiple stars and their relation to star formation. Recently. iRaghavan et al.l (120101) 



confirmed that stellar multiplicity decreases with stellar mass (or later spectral type) in the vicin - 



ity of the Sun, a finding that was already suggested earlier (e.g. ISiegler et al 



2005 



Lada 



2006). 



Because the stellar population in the solar neighborhood is dominated by late-type M-stars, this 
implies that about half of the solar-like main sequence stars in the solar neighborhood occur in 
multiple systems. 

It is likely that the multiplicity rate of stars is establ ished in the early phases of star formation or 



during the dynam ical evolution that occurs afterwards ( [Goodwin et alj2004tlDelgado-Donate et al 



2004 



Bate 2000i. And it is plausible that both initial and environmental conditions constrain mul- 



tiple star formation, and play a vital role in explaining their evolution over time, and their relation 



to the observed MS stellar multiplicity (Mathie u et al 



2008 



Kohler et al. 2008 



2000; 



Sterzik et al 



2003 



Connellev et al 



Kaczmarek et alj|201 lb . It is therefore essential to probe the stellar mul- 



tiplicity as early as possible, to gain a full understanding of the formation and the dynamical evo- 
lutio n of multiple sys t ems. 



tltipl. 

3Ji 



Ghez et all (1993): 



Leinert et 



Duchene et al 



(2004 



al 



Simon et al 



(1995); 



Patience et al 



(2002); 



Beck et al 



119931) : 

2007) all noted differences in the BF of pre-main-sequence (PMS) 



(|2003j)and 

stars in different SFRs. In particular, the multiplicity of T Tauri stars with ages < 10 6 yrs, appears 
to be higher by a factor of 2 in less dense regions such as Taurus comp ared to dens er SFR such as 
Ophiuchus (e.g. 



2002; 



Simon et al. 1995 



Duchene et 



allbooi . 



Beck et al. 2003 



Ratzka et al. 2005) or the field (Patience et al 



Towar d even younger ages , radio continuum (e . 



and near- and mid- Infrared (IR) (Ba rsonv et al 



2005 



Haisch et al 



2002, 



Reipur th et al 



2004 



. 2004) 



2006) investiga- 



tions found additional evidence for an even higher multiplicity among Class I and Flat-spectrum 
(~ 2 - 5 x 10 5 yrs) protostars. 

On the other hand. iMaurv et al.l (|2010) failed to detect companions in a small sample of 5 pro- 
tostars at separations 50 <a <5000 AU (a being the semimajor axis) in a millimeter study of 
se lf-embedded, youn g (Class 0) objects (~ 10 4 - 10 5 yrs). Taken together with the observations 



of Loonev et al 



(2000), the authors argue that multiplicity in the separation ranges from 100 to 600 
AU would only be defined in a later stage of star formation (namely after the Class phase) and 
that early multiplicity in pristine systems may not be as ubiquitous and primeval. 

However, all multiplicity studies in the earliest phases of stellar evolution are plagued with 
relatively small number statistics. In particular, very little is known about companions of embed- 
ded protostars at the sub-AU separation scales, a regime impossible to address with conventional 
imaging techniques. Spectroscopic Binaries (SB) may provide add ition al important constrain ts on 



the star-formation mechanism itself. 



Tokovinin & Sm ekhov (2002) and 



Tokovinin et al 



(2006J) es- 



timated the SB frequency in MS stars and inferred that 65% of their sample of 165 spectroscopic 
binaries were members of higher order multiple systems, often found in a hierarchical configura- 
tion. Remarkably, the frequency of SB in multiple systems strongly depends on the orbital period 
of the SB: 96 of the SBs with orbital periods shorther than 3 days are in multiple systems, while 
this rate drops to 30% for SBs with orbital periods longer than 13 days. The higher frequency of 
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spectroscopic system within triple or higher order systems suggests an imprint of the formation 
mechanism itself. Spectroscopic pairs may have lost their orbital angular momentum owing to the 



( Kozai 


1962; 


Sterzik et al. 


2005: 



Fabrvckv & Tremaine 



2007). This interaction has the effect of 



tightening the spectroscopic binary orbit, while the tertiary would eventually be evacuated to an 
outer region of the system. The SB may therefore hold an important fossil record of (prevailing) 
initial conditions during their formation. However, this scenario is not conclusive by itself, and it 
is not clear at which evolutionary stage this process operates efficiently. Observations are required 
to better constrain the detailed mechanisms of how tight binaries are form e d at v er y young ages. 



(1986); 



The first successful PMS SB detections we re performed by iHerbig 



Mathieu et al 



(1989) and 



(1977); 



Hartmann et al 



Mathieul 11994) in the Taurus and Ophiuchus SFRs with a radial 



velocity precision limited to ~lkms _1 by optical high-resolution spectroscopy. Since then, improv- 
ing instrumental capabilities and b etter radial velocity precision allowed to address PMS sources 
and the lower mass regim es (e.g. iMelol l2003t ICovev et alj|2006t iKurosawa et al.ll2006t IJoergens 



2006, 



2008 



Prato 



2007). Only recently, infrared (IR) high-resolution spectroscopy allowed the 



search for SB among the youngest, embedded PMS and the very low mass stars (e. g. 



Blake et al 



2007 



dFigueira et al 



2010; 



2010b 



Figueiraetal. 



Blake et al. 



2010b 



201 



0; 



Rice et al 



Crockett et al 



2010; 



Prato 



2007 



Crockett et al. 



201 ll) . In particular 



201 ll) reached a precision of less than 50 



ms in PMS stars and very low mass stars. Even the planetary mass regime may be in reach using 
high-stability, high-precision IR spectroscopy. 

Our own work aims to find SBs in a sample of 38 Class 1/ II protostars in the SFR p Ophiuchus. 
We performed a high-resolution IR spectroscopic survey to derive p recision radial velocities (RV) 
using the NIR spectrograph CRIRES at the VLT ( Kaeufl et al]l2004 ). We describe our sample, the 
data analysis methodology and results in Sect. 2. Sect. 3 discusses the RV derived in the context of 
the p Ophiuchus SFR and its implications on the formation of SB at very young ages in general. In 
Sect. 4. we draw our conclusions. 



2. Observational aspects and analysis 

2.1. Observations and data reduction 



The arrival of CRIRES (Kaeufl et al. 



2004), a high-resolution near- and mid-IR spectrograph 
mounted on Antu at the Paranal Observatory allowed for the first time to acquire radial veloc- 
ity measurement precision of a few ms 1 level in this wavelength window. CRIRES has opened the 
possibility to look for (low-mass) spectroscopic companions around young stars that are otherwise 
too faint to be observed in the optical because of the extinction caused by circumstellar material. 
It also entails improved abilit y to detect dimme r companions to optical bright primaries through 
more favorable flux ratios (e.g. Mace et al]l2009 ). 

This study intends to take advantage of CRIRES to measure radial velocities in a sample of 
Class I/II protostars whose multiplicity at larger scales has already been determined by imaging 
techniques. The main purpose consists in detecting radial velocity variations caused by potential 
close companions and to compare the overall multiplicity frequency in embedded sources as a 
function of separation with those observed for T-Tauri and MS stars from the sub-AU to a thousand 
AU scale. 
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extracted nodded spectrum (A) 




2. 296x10' 2.298x10'' 2.3x 10* 2.302x 10 4 2.304x10* 

A(Angstroms) 



Fig. 1. Process of telluric signature removal in DET4. The extracted nodded spectra of GSS26 (A) 
is divided by the spectrum of the telluric standard stars with the same exposure time and observed 
closest in time (same airmass) (B). Lower panel exhibits 'cleaned' nodded spectrum (C). 



To this aim, a sample of 3 8 Class 1/ Flat spectrum ob 



Infrared multiplicity studies o flHaisch et al.l (12002 , 



2004 



ects was selected from the near- and mid- 



2006). The main selection criteria was the 



apparent /T-band magnitude of the targets. As a compromise between the telescope time required 
to reach the desired S/R and the size of the sample, this limit was set at K < 1 1 . In addition, we 
selected proto-stars from 3 different SFRs to allow us a preliminary investigation of the role of the 
environment in our measured spectroscopic binary fraction. 

In CRIPvES the incoming light enters the instrument and is first directed into a pre-disperser, 
after which it is dispersed in a 31.6 lines/mm echelle grating. It delivers high-resolution spectra 
(up to ~100 000 with the a 0.2"-slit) in the 960 - 5200 nm interval in wavelengths windows of 
~ A/70. The observed light is finally projected into four Aladdin III detectors (4096 x 512 pixel) 
separated by ~ 250 pix. The combination of the size of the detectors and the high resolution of 
the instrument limits the wavelength coverage of CRIRES, which is considerably smaller than 
that of other spectrographs. Therefore the first task was to find out the best set-up among the 200 
available to obtain the most accurate radial velocities as a function of the signal-to-noise ratio. 

In order to estimate the S/R necessary to reach the required precision, we investigated the 
amount of information (i.e., number of available spectral lines and their intensity) by simulating 
all CRIRES set-ups as described in iFigueira et al.l (1201 Obi) . In a nutshell, the estimated error on 



the radial velocity is given by e(V r ) °c R [(S/R) < I > n^ 2 ] -1 , where R is the instrumental 
resolution, e is the photon noise error, S/R is the signal-to-noise ratio, < / > is the average depth of 
absorption lines, and n\ is the number of lines entering in the computation of the CCF. In order to 
simulate all CRIRES set-ups, we used t he IR solar spectru m obtained with the Fourier Transform 
Spectrograph provided by the NOAO (Wal lace et al.lll993l) . This solar spectrum was convolved 
with a Gaussian function and trimmed t o match the resolutio n and wavelength coverage of the 
various CRIRES set-ups (for details see iFigueira et al.ll2010bl) . Assuming that only stellar lines 
deeper than a certain threshold contribute to the CCF, the above expression was applied to all 



4 



P. Viana Almeida et al.: Finding proto-spectroscopic binaries 

CRIRES simulated spectra and the quantity e(Vr)/S/R was derived as a function of the set-up. 
This ratio allowed us to compare the radial velocity error (in a relative way) among the set-ups and 
to choose the most s uitable one. A few set-ups were identified in the H- and K-band (see Fig. 1 of 



Figueira et al.l l2010b). Because of the unknown degree of obscuration of our sample, we selected 
the K-band set-up. 

To derive the absolute error and therefore be able to transform the ratio e{Vr)/S/R into 
exposure times, we had to know the absolute error on V r as a function of the S/R for one set-up. 
Our exposures times were then computed assuming that an S/R of 20-30 in the setting chosen 
here would give a 150m/s error for the set-up former 24/-l/i setting. In practice, this 
estimate turned out to be too optimistic. As a consequence, exposure times for several targets were 
underestimated. Accordingly, we obtained 15 protostars with robust RV measurements from the 
original 38 stars proposed and only 7 sources with multiple RV measurements at sufficient S/R. 
This paper focus on these 7 multi-epoch targets from the p Ophiuchus SFR, for which we detected 
photospheric absorption features. In Table Q] we present the properties of the full sample. An 
assessment of the evolutionary status of each sources is presented in the column SED of Table Q] 
It follows the usual YSO classification based on the change of the slope of the SED as a function 
of age. Rising SEDs with positive slopes toward mid-IR usually depict the presence of protostellar 
envelopes that are generally more massive than the protostar. These protostellar envelopes tend 
to disappear in the subsequent phases of evolution, namely the Class I phase. As a consequence, 
the SED slope flatte ns (flat-spectrum sou rces) and then becomes negative towards the Class II or 
T-Tauri phase, (e.g. Motte & AndrelboOl ). 



Our sample was observed between April 2008 and February 2009 with a set-up covering a 
wavelength range of 2254.2 to 2304.7 nm. The center wavelength reference in detector3 (DET3) 
was /J2287.1 nm. The CO (2-0) bandhead at 2293. 5nm was registered on detector 4 (DET4). 
Owing to the lack of reference stars close to our science targets, our observations were carried 
out in the seeing-limited mode. Observations were collected under good seeing (~ 0.8") in 2 
nodding cycles (i.e., 4 total independent exposures)^, with a 0.3" wide slit allowing a resolution 
of ~ 60,000 which we found to be the best compromise between the amount of flux entering the 
spectrograph and the ability to resolve the lines of the object and those from the atmosphere. Final 
S/N ratios given in TableQ]were computed on small featureless regions. 

The approximate range of masses of our 7 protostellar targets was taken from 



Natta et al 



(2006) and assumed to be within the interval from 0.18 M Q to 1.4 M Q . Assuming that semi- 
amplitudes of 500m/s can be measured at a 3cr level, we are able to detect companions ranging 
from several jupiter masses with periods of one or a few weeks to equal-mass companions (from 
brown dwarfs to solar-type stars) with periods of several y ears. A more detailed study of the 



detectability rate similar to the one performed bv lMelol (120031) will be carried out at the end of the 
survey. 

In addition to the protostellar data we also collected spectra of telluric standard stars (fea- 
tureless early-type B stars) observed with the same airmass and instrumental setup as our science 
targets. Multi-epoch spectra of 3 radial velocity standard stars were collected during the timespan 



1 In one nodding cycle two exposures are taken in total, one at position A and another at position B. In the 
second nod cycle, a third exposure is taken again at position B, then the telescope is moved back to position 
A where the 4th exposure is taken. 
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Fig. 2. Extracted spectra of our Class-I stars as observed in Detector 4. The data shown here were 
smoothed for clarity, continuum-normalized and Doppler-corrected to a null radial velocity. 



of our observations in order to study the precision of our RV measurements. 

Data were reduced using an optimized IRAF-based pipeline to ensure a uniform and homoge- 
neous reduction for all targets in the different epochs. The reduction steps in our pipeline were: 1) 
correction for CRIRES nonlinearity using ESO on-line coefficients and construction of bad-pixel 
masks, 2) dark-current subtraction from the fiatfield images and collapsing into a master flat, 
3) division of science spectra by these final fiatfield images, 4) su btraction of op posing nodded 
spectra for dark-current removal, 5) optimal extraction using the iHornel d 19860 algorithm, 6) 
telluric signature removal by dividing each extracted spectrum by the spectrum of the telluric STD 
stars with the same exposure time and observed closest in time (same airmass) (Fig.[TJ. 

The position of the CCF peak strongly depends on the average center of all stellar lines used 
to computed the CCF. Because telluric removal could have changed the final line shapes, it is 
essential to find the most similar and adequate telluric spectra so as not to change the protostellar 
line shapes fundamentally. A test was carried out to assess the impact of the use of different 
observed telluric spectra in the derived RVs. First, we performed for a specific date telluric line 
removal with the observed telluric standard and the combined science spectra of protostars where 
no stellar signal was detected (i.e., only the atmospheric lines are imprinted in the data). The 
main conclusion was that small fluctuations (< 50ms~') on the final RV could indeed be detected. 
These discrepancies were found in tellurics that were observed >2 hours after the science target. 
We found no relevant differences in the final RV when the telluric removal was performed with 
tellurics observed within that time difference. We conclude that telluric removal as explained here 
has no major impact on the final RV. 

The extracted spectra were wavelength-calibrated using telluric absorption line s. The use of 
tellur i c features to build a wavelength solution i s now a well-documented technic 



1985 



Blake et al 



2007 



Figueira et al 



2010b 



Blake et al. 



2010; 



Crockett et al. 



ue (Caccin et al 



2011) and it is 
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highly convenient when gas cells are not av ailable as was the ca se of CRIRES a t the time of 

J20ld) for an 



the observations. We refer the reader to (e.g. IFigueira et al 



2010a) or 



Blake et al. 



2010aJ), in 



extensive explanation of the robustness of this technique. The study of (Fig ueira et al. 
particular, showed that 02 atmospheric lines can be stable down to 10 m/s over long periods of 
time (6 years). 

Basically, wavelength solution construction consists in finding and fitting (in each observed 
spectrum) the center of telluric lines with low-order polynomials to match their theoretical 
value as given by online databases. In our particular setting, we took advantage of the wealth 
of CH4 atmospheric telluric lines to build suitable wavelength solutions for each extracted 
spectrum. This set of CH4 lines has already been used in previous s tudies and was shown t o 



deliver precise RVs in very low-mass stars and in brown dwarves (Blak e et al 



2007 



2010). 



The reference wavelength of each spectral line was taken from the HITRAN database with 
a typical accuracy < 50ms _1 , which is well-suited to our purposes, giving the radial velocity 
amplitudes expected by the presence of stellar companions. Spectral fea tures of protostars are 



known to be considera bly broader than MS stars or even class-II stars (iGreene & Lada 



1997 



Doppmann et alj 120051) . because of their typically high vsini. For this reason, we expected to 



find a high level of blending between spectral features of telluric and stellar origin that could 
introduce significant systematic errors when building the wavelength solution or determining 
the RV measurements. To minimize possible systematic errors arising from the overlap of the 
telluric with the stellar signal, we combined all nodded images from each target's observation 
into a single frame for building the wavelength solutions. This procedure considerably improved 
the S/N of the final spectra, which is therefore ensured a better identification and fitting of 
telluric lines and their theoretical centers. This steps assumes that the spectral characteristics 
of the atmosphere did not vary between nodded positions (i.e., in time-scales of <4 min in K-band). 



2.2. Data analysis 



The RV s in this study were derived using a slightly different version of the IFigueira et al.l 
(2010b) p ipeline, a two -dimensional (2D ) cross-correla t ion fun ction (CCF) based on TODCOR 



( Maze h & Zuckei 



1992). The pipeline of IFigueira et al.l d2010b ) was specifically designed to de 



liver the RV of the target relative to the zero-point established by the Earth's atmosphere. 

In a 2D cross-correlation using atmospheric lines as a reference, the RV of the atmosphere 
(RVjfcy) and the RV of the star (RV Jtor ) were obtained from the same spectrum by cross-correlating 
the observed spectrum with a combination o f two CORAVEL-ty pe masks (i.e., the absorption lines 
are represented by box-shaped lines, see iBaranne et al.l ll996). one mask containing the telluric 



lines and another the stellar lines. As mentioned in the previous section, the lines present in the 
atmospheric mask use the HITRAN wavelengths as reference. The stellar masks were built using 
the synthetic spectra of PHOENIX for stellar type G8V, K2V, K5V, Ml and M6. The width of the 



box-s haped line was set by convolving the synthetic spectra with Gray's rotational profiles (Gray 



1992) with the rotational velocities chosen to adjust the line width of the observed spectra. The 
limb darkening coefficient was found to be irrelevant so we assumed a standard value of 0.6. 
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RV standard star: Gl 406 




RV Maskl (star) 



Fig. 3. Resulting 2D CCF for the star G1406 and telluric lines using an M6 stellar mask and an 
atmospheric mask, respectively. The best correlation is found for the atmosphere (RV Mask2) and 
the star (RV Maskl) at the lowest point of the function in the z-axis. 



We stress that different combinations of spectral type (i.e., logg and T e ff), rotational velocity 
and veiling would well fit the observed spectra equally. For this reason, the process above was 
useful to set the line center and width and to generate the box-shaped mask that was used in CCF 
computation. No attempt was made to derive the target's T e ff and v sin i. To test the sensitivity of 
our radial velocities with the stellar mask, we built a new series of stellar masks using the signal- 
to-noise spectrum of our RV standards HD129642 (K2V), G1406 (K5V) and HD 105671 (M6). 
We found that changing the spectrum to build the masks induced RV variations below 100 m s _1 , 
which is within our photon noise. 

The best solutions for the RV SJ t v and RV J/flr were found by fitting the 2D CCF surface (as shown 
in Figures with a two-dimensional Gaussian function. The coordinates of the m inimum give RV 



Baranne et al. 1996) 



and By s tar- (for a description of the ID CCF technique see 

In our s pecific case, the most im portant aspect of our adapted pipeline compared to the original 
' Figueira et al. ( 20 1 Obi) is that it uses two spectra: the one produced after the telluric 



2D CCF of 



removal (e.g. spectrum C in Fig. [TJ when computing the stellar CCF against a stellar mask, and 
the original one containing the stellar and the telluric lines [e.g. spectrum A in the same figure) 
when computing the atmospheric CCF against an atmospheric mask. Although it uses two spectra, 
the process is still performed simultaneously, so that the final product is still a 2D function of the 
velocity shifts of the two templates (see Fig . [3j . All wavelength solutions and derived RVs were cor- 



rected for the solar system center of mass following the references given bv lBretagnon & Francou 



( 1988). Barycentric Julian dates were also retrieved according to the same source. The final stellar 
velocity is then given by RVfi„ a i tSlar - R V star - RV s k x + BERV, BERV being the barycentric Earth's 
radial velocity on the UT date on which each target was observed. 

We note that during telluric removal the final S/N of the cleaned spectrum (C in Fig. [TJ de- 
creased to about ~20. The S/N of individual nodded positions was often low. Therefore, with the 
exceptions of the sources GSS26 and GY23 and the RV standards, we combined the nodded spectra 
to form a single spectrum in order to increase the final S/N. 

In addition, we also included an algorithm (analog to the continuum task from the IRAF pack- 
age) to better fit the noisy continuum around each spectral line used in the cross-correlation. This 
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last procedure had the effect of "cleaning" the final shape of the CCF from spurious spatial fre- 
quencies introduced by the incorrect fitting of the continuum level in the noisiest spectra. 



2.3. Results 

The RV measurements of stars that possess multi-epoch observations are shown in Table[T]and are 
plotted in Fig.[4]and[6] Final RVs come from the 2D CCFs calculated using the spectral information 
in the wavelength region of the CO Av = (0-2) bandhead window sampled in DET4 in the chosen 
CRIRES setting (~ 14 absorption lines). The errors quoted in Table[T]were computed as follows. 
For each independent exposure (nodded spectrum), an RV was derived for each of the 14 lines 
in the region of the CO Ay = (0-2) bandhead. A weighted mean was calculated using the line 
depth as weight. An average deviation of the 14 RVs was then computed. The final cr quoted in 
TableQ]is then the average cr divided by N^ 2 , where N is the number of the total nodded spectra. 
Errors for GSS26, GY23 and the RV standards (HD129642, HD 105671, G1406) were computed as 
described above. For the remaining targets where we had to combine all available exposures, the 
errors presented in Table[T]are the RMS of the RVs computed using the 14 lines. The errors range 
from around 0.015 (in GSS26) to 1.20 kms -1 in the noisier stars (S/N-15). 

We caution t hat the low errors obt aine d for some stars in Ta ble [T] are probably underestimated. 



According to lFigueira et al 



(2010a) and 



Figueira et al.1 yOlOb), the precision one can obtain using 



the telluric lines as zero-point is better than 10m s _1 . This is about what we obtain on the RV 
standards. On the one hand these low uncertainties are telling us that the systematic errors (target 
centering, detector characteristics, etc.) are not dramatically higher than those fund on the time- 
scale of our measurements (100 days or so). On the other, they might simply reflect the low-number 
of measurements obtained for our targets. A more realistic estimation of the errors will therefore 
be obtained later in our survey. 

To assess the precision obtained by the method presented here, we computed the RV for the RV 
standards G1406, HD129642, and HD 105671. For G1406, we obtain 19.730 kms 1 with an internal 
dispersion of 0.010 kms" 1 in all nodded images. Our results is very close to the published result 
provided in lwhite & Basril (120031) of 19. 175 kms" 1 and to the 19.18+0.03 km s" 1 obtained over 11 
measurements (Forveille, private communication). Quite striking is the RV = -4.59 kms -1 and the 
low dispersion (~0.08 kms 1 ) obtained in the four different epochs in th e star HD129642. R esults 
are plotted in Fig [4] An RV precise value of -4.60 kms 1 is given by iGontcharovl (|2006). The 



four epochs obtained for HD 129642 coincide with the temporal window of most of our Class I/II 
sources. That they are able to reproduce HD 129642 RV during such an extended period confirms 
the stability and precision of the results. 



3. Discussion 

3.1. Origin of the CO: photospheric or circumstellar 

Protostellar light is often reprocessed in the circumstellar environment and becomes contaminated 
with emission or absorption from different origins. P articularly in the near-IR the circumstellar 
disk is thought to play a central role ( Calvet et al.l 199 1 ) in altering the stellar signature because the 
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S -4.6 - 



2.4546x10 s 2.45465x10 s 2.4547x10" 

HJD(days) 



Fig. 4. RV measurements for the STD star HD129642. Error bars come from the relation «T vra d>* 
N~ 1/2 being cr vrad the dispersion of the RV determined in each nodded image and N the number of 
individual measurements. 



temperatures of the excited material emit in the same wavelength as the protostellar pho tosphere. 



Fo r the CO regions, for instance, dedicated stu dies in TTS with accretion disks (e.g. 



1993 



Naiita et al 



Carr et al 



1996; 



Luhman & Rieke 



1999) have uncovered observations of CO overtone 



bandheads. The presence of inclined circumstellar disks can affect the shape of observed CO spec- 
tral lines that can exhibit anomalous double peaking profiles or depict a wider velocity field gradi- 
ent. These features can be o bserved mainly in emission, but can also be found in absorption (see 
e.g. iHoffmeister et al.ll2006[ and references therei n) and are apparently associated with the most 
active PMS sources such as the Fu Orionis objects (IHartmann et al CO absorption lines can 

be produced in a circumstellar disk when its t emperature profile depicts a negative vertical gradient 
very much like a stellar photosphere (see e.g. 



Kaufl et al. 2005). 



Nevertheless, because the dynami cs of the material in circumstellar disks present Keplerian 
velocity profiles (ICasali & Eir oa 1996), the CO absorption of circumstellar nature drastically con- 
trasts with the CO absorption profiles of photospheric origin. High-resolution observations are 
usually able to separate both components. 

In Fig. |5] we show that our protostellar spectra can be reproduced by rotationally broadened 
photospheric synthetic spectra. As an example in this figure, we rotationally broadened 5 solar 
metallicity PHOENIX stellar models of Spectral type G8V, K2V, K5V, Ml and M6 (T eff ranging 
from 5650K to 3050K) with the same resolution as our observations and overplotted them with the 
spectrum of the source GY23. The match between the synthetic and the observed profiles is evi- 
dent. Using different sets of models we were able to extend this result to all 7 protostars presented 
in this work. For this reason we believe that the lines analyzed here are generated in a photospheric 
environment and that the RVs found are of stellar origin. 
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2.297 2.298 2.299 

fj.ni 



2.297 2.298 2.299 

j.im 



Fig. 5. Observed protostellar spectra that we were able to reproduce with rotationally broadened 
stellar photospheric profiles. Right panel: Overplotted are some synthetic spectra (with a fixed T e ff 
of 3649K) (solid lines) with different vsini rotations with the spectrum of the protostar GY23 on 
one of the dates it was observed. Left panel: The same, but with a fixed vsini of 17kms~' and 
different temperatures. 



3.2. Radial velocities 

From the RV measurements presented in Table Q] we can see that the range of values obtained 
(from ~ -5.5 to -17 kms~') seems to indicate that the sources studied here are indeed from the p 
Ophiuchus cloud. All targets, with the exception of VS SG18, are within 3 cr of the chara cteristic 
velocities exhibited by pOph sources, -6.3 + 1.0 kms" 1 (IKurosawa et al 



Prato 2007). 



The phenomenon of accretion and stellar activity in protostellar ages has the effect of altering 
spectral line-shapes. These non-periodical variations change the l ine-shape and induc e RV varia- 
tions that can reach amplitudes of the order on the -2-3 kms -1 {e.g. 



Alencar et al 



2005). Therefore, 



although it is possible that some of the observed RV variations might be caused by stellar activity 
and accretion, the amplitude of the variations (^ 4-6 kms 1 ) observed for VSSG18, GY23, and 
GY5 1 suggests that those variations could have been caused by a spectroscopic companion. 

Moreover, for the other 3 proto-stars observed with lower S/R (IRS34, L1689SN02 and 
WL17) low-amplitude RV variations were detected, suggesting that the RV variations found among 
VSSG18, GY23, and GY51 might well be produced by companions. Because the effects of this 
accretion-driven RV modulation do not persist for long (see ibid), the only way to confirm the 
nature of the observed variations is to assess the eventual variations of the orbital parameters. In 
particular, variations of the semi-amplitude and/or the phase w ould indicate that the RV modula- 
tions are likely to be accretion-driven (e.g. Huelamo et alj|2008 ). 



Finally, concerning VSSG18 alone, al l RV measurements are well above the 3 cr of the mean 
of the velocity distribution of pOph from IKurosawa et ajj (120061) . VSSG18 average velocity is of 



-15.873 kms l , a clear but not extreme outlier from the local cloud. Based on its spectral energy 
distribution and spectral index from previous imaging studies it seems that VSSG18 is indeed a 
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Table 1. Barycentric Julian dates of each target, derived RVs, and final dispersion for each date it 
was observed. In the upper table the results from each nodded image are presented. In the lower 
panel we show the RV measurements in images of sources that yielded an insufficiently S/N after 
telluric removal and that had to be combined for better precision. 



Source 


MJD 


j n —In 

vrad [km s '] 


<cr> 


Source 


MJD 


j n —In 

vrad [km s '] 


<cr> 


nodded images 
















GSS 26 


2454584.63565 
2454584.63909 

r\ A C A C Ct A f "\ — 1 A C 

2454584.63745 


-6.97198 
-6.93852 
-6.94756 




HD 129642 


/-> ,( r* A r" t \ A if f\ l—f o a 

2454584.62788 
2454584.63176 

r\ a c a c c\ a f i "i n — j 

2454584.62387 


A A d A f\ i*\ 

-4.48490 

A c f\r\ a f\ 

-4.59242 
-4.22596 






r\ A C A C Ct A if A f\Ct — 1 

2454584.64087 


-6.96722 


0.01 




f\ A C A C C\ A f 1 t'\ — 1 if 

2454584.61976 


,( C f\f\ A 

-4.59242 






2454595.81134 


-6.93863 






2454584.63573 


-4.60915 






2454595.81313 


-7.03562 






2454584.63931 


-4.26827 


0.09 




2454595.81492 


-6.95636 






2454595.76527 


-4.48264 






2454595.81664 


-6.89839 


0.03 




2454595.76193 


-4.71534 






2454700.50155 


-6.98628 






2454595.75729 


-4.67921 






2454700.50334 


-6.99655 






2454595.75303 


-4.64180 






2454700.50501 


-6.94632 






2454595.76937 


-4.72930 






2454700.50681 


-6.86636 


0.03 




2454595.77303 


-4.82873 


0.06 


GY23 


2454584.66269 
2454584.66334 
2454584.66383 


-11.74129 
-11.91061 
-11.61631 
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r\ a c a if r\ if r~ i i / \ t i 
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-4.59595 
-4.49518 
-4.63937 
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-1 -1 f\f\Clf\if 

-11.92826 
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f\ A C A if r\ if c t \ — J -i f\ 

2454696.50713 


-4.34328 






2454595.80848 


-6.86636 






f\ a c a if r\ if c f\ f\ -\ r\ 

2454696.52319 


-4.60013 






f\ a c A c r\ r~ 1 1 i\ — 1 1 1 f 

2454595.80786 


-6.99655 






f\ A C A if r\ if c f \ — 1 f\ f\ 

2454696.52733 


-4.75527 


0.07 




2454595.80737 


-6.94632 






f\ A C A — I 1 <-\ C* r\f\ f\ -\ 

2454712.50231 


-4.80776 






/4j4jyj.oUO 12. 


a com 1 
-O.oZUl 1 


f\ f\A 




1 A <AH A 1 <f\AQf\ 
Z4j4 / 1Z.JU48U 


A CQO/;'! 

-4.jyzoj 






2454710.51703 


-10.14933 






2454712.50349 


-4.68273 






2454710.51638 


-9.97590 






2454712.50343 


-4.72138 






2454710.51590 


-10.23329 






2454712.50391 


-4.44132 






2454710.51526 


-10.11078 


0.05 




2454712.50429 


-4.62820 


0.07 


G1406 


2454890.64318 


19.73374 




HD 105671 


2454795.13436 


-4.20193 






2454890.64215 


19.72979 






2454795.13447 


-4.22028 


0.01 




2454890.64125 


19.74171 














2454890.64024 


19.70016 


0.01 










combined images 


GY51 


2454584.67379 


-10.75198 


0.22 


L 1689 SN02 


2454699.66023 


-7.76361 


0.76 




2454595.83890 


-9.44251 


0.35 




2454873.87196 


-7.26417 


0.56 




2454597.72588 


-6.26190 


0.25 




2454884.87905 


-5.76223 


0.92 




2454710.527356 


-10.18371 


0.73 










VSSG 18 


2454597.88737 


-15.89123 


0.28 


IRS 34 


2454597.80272 


-7.36514 


0.63 




2454612.72451 


-13.73072 


0.32 




2454612.80033 


-7.17129 


0.80 




2454733.55090 


-17.12639 


0.44 




2454712.56017 


-7.25198 


1.22 


WL 17 


2454597.75984 
2454610.76480 


-5.62187 
-6.12428 


0.65 
1.10 











proto-star. Its deviant RV could suggest that this is long -period SB1 system. 



3.3. Multiplicity and star formation scenarios 

The use of CRIRES along with the technique described in the present work allowed us for the first 
time to carry out a systematic search for spectroscopic companions at early stages of their formation 
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Fig. 6. Radial velocity measurements of the 7 protostars of our sample. Error bars in GSS26 and 
GY23 are calculated from the relation «r vra d>* N -1 ^ 2 being cr wa d the dispersion of the RV deter- 
mined in each nodded image and N the number of individual measurements. In all others targets 
error bars were computed as explained in Sect. 2.4.. 



(below about lMyr or so). Although preliminary and covering a small number of objects (7 objects 
only), the present paper brings already some interesting aspects concerning the configuration where 
those systems are found. 

Out of 7 sour ces, 4 ar e known to be in multiple visual systems, n amely, L168 9NQ2 



(Haisch et al 



1978 



Haisch et al 



2004), IRS34 (Bar sonv et al 



2002 



2005), GY51 (lHaisch et al. 



2006), and GY23 (Elias 



2004), whereas the other 3 targets (GSS 26, WL17 and VSSG18) 



did not have any detected companions in previous imaging surveys (lHaisch et al 



2004 



2006 



Barsonv et al. 2005 



Duchene et al. 2004). 



GY5 1 and GY23 resemble some young multiple s ystems found in the literature. The former, is 



maybe a non-hierarchical system such as HD 34700 (iSterzik et al 



2005) whereas the latter looks 



like to be in the typical hierarchical configuration with a tight inner binary wi th an additional 



comp anion which probably caused the inner orbital evolution of the inner pair (ITokovinin et al 
20061) . 



The detection of a potential spectroscopic VSSG18 system without visual companion detected 
is interesting. If we believe that one possible (or maybe the only possible) way to form spectro- 
scopic pairs is through the existence of an additional companion, systems such as VSSG18 are 
valuable because, if confirmed, they might set time constraints for orbital evolution to occur. Also, 
they could be interesting targets for AO-assisted campaigns to look f or more close c ompan i ons 



When we com binethe results of our stud y with those of 
J2005L and 



Haisch et al 



Barso nv et al. 



Duchene et al 



(2004 



2006), 



(2004), there were 2 single systems (S) amongst the 
YSO analyzed systems, 3 binaries (B), 1 triple (T) (1 spectroscopic binary + 1 visual companion), 
and one quadruple (Q) (1 spectroscopic binary + 3 visual companions). The binary fraction (BF) 
of this small sample of protostars (B + T + Q)/ (S + B + T + Q) is ^71%. This preliminary result, 
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although based on such a small sample, agrees surprisingly well with the notion that multiplicity 
is very high at young ages and therefore it might be a product of star formation. We advise caution 
in the interpretation of this final BF, however, because it is likely affected by low number statistics 
and it requires confirmation with a larger sample of protostars. 

4. Conclusions 

We have performed a spectroscopic multi-epoch survey of 7 known protostars in the pOphiuchus 
star-forming region in the near-IR. Our main goal was to derive the close-binary frequency of 
embedded sources where multiplicity is already known from imaging techniques. We successfully 
derived radial velocities of Class I/II sources with an unprecedented precision in the spectroscopic 
study of protostars. Depending on the orbital period, our method allows for the detection of 
companions within the planetary-mass domain around embedded sources. 

We find tentative evidence for the existence of three spectroscopic multiple systems out of 
the total 7 protostars analyzed. With an internal precision in the range of 0.02 to 1.22 kms -1 
we detected clear RV variations on the order of 4 to 6 kms 1 . When we combined our results 
with those of other studies, we obtain a binary fraction of ^71.5% for these systems, which is in 
line with the idea that multiplicity is high at protostellar ages. The existence of such young SBs 
strengthens the notion that dynamical evolution has already taken place in the Class 1/ II phase 
(ages ~ 10 5 yrs). 

Future observations are needed to enhance the sample statistics and will aim to derive the 
orbital solutions for these SBs if their multiple system status is confirmed. 
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Table 1. From left to right we present the star-forming region of each target, the coordinates, observation 
date, integration time of our observations, K-magnitude, the number of nodding cycles observed, the S/N*, 
the signal to noise ratio of the extracted spectra before telluric removal, the spectral energy distribution (SED), 
the reference from where the SEDs were collected and an alternative name for each source, if available. 
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Class 


Hai04 , Hai06 


V* FK Cha 










2008-12-27 


30.0 




2 


30 
















2008-12-28 


30.0 




2 


20 
















2009-02-23 


30.0 




2 


20 








ISO-Cha I 26 


Cha 


1 1 08 04.00 


-77 38 42.0 


2008-04-28 


45.0 


8.2 


2 


30 


Class I 


Hai04 , Hai06 


HD 97048 2 










2008-12-28 


45.0 




2 


25 
















2009-02-23 


45.0 




2 


30 








Cha I T32 


Cha 


11 08 04.61 


-77 39 16.9 


2008-04-28 


30.0 


6.1 


2 


25 


Class 


Hai04 , Hai06 












2008-12-28 


30.0 




2 


25 
















2009-02-23 


30.0 




2 


25 








Cha I T33B 


Cha 


11 08 15.69 


-77 33 47.1 


2008-04-28 


30.0 


6.9 


2 


30 


FS 


Hai06 


Glass la 










2008-12-28 


30.0 




2 


35 
















2009-02-23 


30.0 




2 


30 








Cha I C9-2 


Cha 


11 08 37.37 


-77 43 53.5 


2008-04-28 


45.0 


8.6 


2 


20 


Class 


Hai04 , Hai06 












2008-12-28 


45.0 




2 


20 
















2009-02-10 


45.0 




2 


25 
















2009-02-26 


45.0 




2 


20 








ChaICl-6 


Cha 


11 09 23.30 


-76 34 36.2 


2008-04-28 


45.0 


8.4 


2 


25 


Class I 


Hai04 , Hai06 


CCE98 1-76 










2008-12-28 


45.0 




2 


30 
















2009-02-26 


45.0 




2 


25 








Cha I T41 


Cha 


1 1 09 50.39 


-76 36 47.6 


2008-04-28 


30.0 


7.0 


2 


30 
















2008-12-28 


30.0 




2 


25 
















2009-02-26 


30.0 




2 


25 








Cha I T42 


Cha 


11 09 54.66 


-76 34 23.7 


2008-04-28 


30.0 


7.0 


2 


20 


Class 


Hai04 , Hai06 


V* FM Cha 










2009-02-10 


30.0 




2 


20 
















2009-02-26 


30.0 




2 


20 








Cha I T44 


Cha 


11 10 01.35 


-76 34 55.8 


2008-04-28 


30.0 


6.4 


2 


20 


Class 


Hai04 , Hai06 


CHXR 44 










2009-02-10 


30.0 




2 


30 
















2009-02-26 


30.0 




2 


25 








Cha I T47 


Cha 


11 10 50.78 


-77 17 50.6 


2008-04-28 


45.0 


8.8 


2 


20 






V* FO Cha 
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Source 


Region 


ff(2000) 


5(2000) 


Date 


DIT (sec) 


Kt map] 


nod.cy 


S/N* 


SED 


reference 


alias 










2009-02-11 


45.0 




2 


30 
















2009-02-26 


45.0 




2 


20 








Cha TT 8 

V^llLL XX (J 


Cha 


12 53 42 88 


-77 1 5 05 7 


2008-05-11 


45.0 


8.8 


2 


25 
















2008-08-17 


45.0 




2 


20 
















2008-08-21 


45.0 




2 


30 








GSS 26 


pOph 


1 6 26 1 28 

1U ± VJ . o 


-24 20 56 6 

Z.T 1 Z.W JU.U 


2008-04-28 


120.0 


9.4 


2 


25 


Class I 


Bar05 












2008-05-09 


120.0 




2 


35 
















2008-08-21 


120.0 




2 


30 








GSS/IRS 1 


pOph 


16 26 21.50 


-24 23 07 


2008-04-28 


120.0 


9.0 


2 


30 


Class I 


Bar05 












2008-05-09 


120.0 




2 


20 
















2008-09-01 

yj y yj i 


120.0 




2 


30 








GY 23 


pOph 


16 26 24.00 


-24 24 49.9 


2008-04-28 


30.0 


7.4 


2 


30 


Class TT 


Bar05 


S2 










2008-05-09 


30.0 




2 


25 
















2008-09-01 

yj y yj i 


30.0 




2 


30 








GY 51 


pOph 


1 6 26 30 49 


-24 22 59.0 


2008-04-28 


180.0 


10.2 


2 


30 


FS 


Hai06 


VSSG27 










2008-05-09 

\jvhj yj^j yj y 


180.0 




2 


30 
















2008-05- 1 1 


180.0 




2 


35 
















2008-09-01 

^uuu yj y yj i 


180.0 




2 


25 








WL 12 


pOph 


1 6 26 44 30 


-24 34 47.5 


2008-04-28 


180.0 


10.4 


2 


15 


Class I 


Bar05 


GY111 










2008-05-09 


180.0 




2 


20 
















2008-09-02 


180.0 




2 


20 








WL 1 S 


pOph 


1 6 27 04 13 


-24 28 30 7 


2008-05-09 


1 80.0 


10.8 


2 


25 


Class TT 


Bar05 












2008-05-13 


1 80.0 




2 


20 








WL 17 


pOph 


1 6 27 06 79 


-24 38 14.6 


2008-05-11 


1 80.0 


10.3 


2 


25 


Class I 


Bar05 


GY205 










2008-05-24 


180.0 




2 


20 








Flias 29 


pOph 


1 6 27 09 43 

±\j — 1 yj y .^^j 


-24 37 1 8 5 


2008-05-1 1 


30.0 


7.5 


2 


25 


Class 


Hai04 Hai06 


GY214 










2008-05-26 


30.0 




2 


20 
















2008-09-03 


30.0 




2 


20 








GY 224 


pOph 


1 6 27 1 1 17 

I VI _ / ±1.1/ 


-24 40 46 7 


2008-05-1 1 


180.0 


10.8 


2 


20 


FS 


Bar05 


WLY 1-43 










2008-05-24 


180.0 




2 


30 
















2008-09-03 


180.0 




2 


25 








IRS 34 


pOph 


1 6 27 1 5 48 


-24 26 40.6 


2008-05- 1 1 


180.0 


10.3 


2 


25 


FS 


Bar05 


GY239 










2008-05-26 


180.0 




2 


25 
















2008-09-03 


180.0 




2 


30 








IRS 37 


pOph 


16 27 17.54 


-24 28 56.5 


2008-05-11 


180.0 


10.9 


2 


25 


Class I 


Hai04 


GY244 










2008-05-26 


180.0 




2 


25 
















2008-09-24 


180.0 




2 


25 








IRS 42 


pOph 


16 27 21.45 


-24 41 42.8 


2008-05-11 


45.0 


8.6 


2 


25 


FS 


Bar05 


GY252 
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Source 


Region 


ff(2000) 


(5(2000) 


Date 


DIT (sec) 


Kt map] 


nod.cy 


S/N* 


SED 


reference 


alias 










2008-05-24 


45.0 




2 


30 
















2008-09-24 


45.0 




2 


30 








WL 6 


pOph 


16 27 21 83 

IV'—' z- i . u *j 


-24 29 53.2 


2008-05-11 


180.0 


10.8 


2 


35 


Class I 


Bar05 


GY254 










2008-09-24 


1 80.0 




2 


30 
















2008-09-25 


180.0 




2 


30 








IRS 43 


pOph 


1 6 27 26 90 


-24 40 51 5 


2008-05-1 1 


120.0 


9.4 


2 


25 




Bar05 


GY265 










2008-05-24 


120.0 




2 


20 
















2008-09-03 

\JVHJ yj y yj^t 


120.0 




2 


25 








IRS 44 


pOph 


1 6 27 28 00 


-24 39 34 3 


2008-05- 1 1 


120.0 


9.7 


2 


30 


Class I 


Bar05 


GY269 










2008-05-24 


120.0 




2 


20 
















2008-09-03 


120.0 




2 


30 








VSSG 1 8 


pOph 


16 27 28.44 


-24 27 21.9 


2008-05-1 1 


120.0 


9.2 


2 


30 


FS 


Bar05 












2008-05-26 


120.0 




2 


30 
















2008-09-24 


120.0 




2 


30 








IRS 46 


pOph 


1 6 27 29 70 


-24 39 16 


2008-05- 1 1 


1 80.0 


10.6 


2 


20 




Bar05 


GY274 










2008-05-26 


180.0 




2 


20 
















2009-02-23 


180.0 




2 


30 








GY 279 


pOph 


1 6 27 30 1 8 


-24 27 44.3 


2008-05- 1 1 


45.0 


9.0 


2 


30 


FS 


Bar05 


IRS47 










2008-05-26 


45.0 




2 


20 
















2008-09-25 


45.0 




2 


25 








IRS 48 


pOph 


1 6 27 37 20 


-24 30 34 


2008-05-11 


30.0 


7.7 


2 


20 


FS 


Hai06 


GY304 










2008-05-26 


30.0 




2 


25 
















2009-02-22 


30.0 




2 


25 








TRS 5 1 S 


pOph 


1 6 27 39 84 


-24 43 16.1 


2008-08-21 


45.0 


8.7 


2 


20 


FS 


Hai06 


GY315 










2008-09-03 


45.0 




2 


25 
















2009-02-22 


45.0 




2 


30 








IRS 54 


pOph 


1 6 27 5 1 70 


-24 31 46 


2008-08-21 


1 80.0 


10.2 


2 


20 


Class I 


Hai06 


GY378 










2009-02-12 


180.0 




2 


25 
















2009-02-22 


1 80.0 




2 


30 








IRS 63 


pOph 


16 31 35 53 


-24 01 28.3 


2008-08-21 


120.0 


9.3 


2 


25 


Class I 


Hai02 Hai06 


WI Y 2-63 










2009-02-11 


30.0 




2 


30 
















2009-02-23 


30.0 




2 


20 








L1689SN02N 


pOph 


16 31 52.13 


-24 56 15.2 


2008-08-21 


45.0 


8.3 


2 


30 


FS 


Hai06 












2009-02-11 


45.0 




2 


20 
















2009-02-22 


45.0 




2 


25 








IRS 67 


pOph 


16 3201.00 


-24 56 44.0 


2008-08-25 


180.0 


10.3 


2 


20 


Class I 


Hai02 , Hai06 












2009-02-11 


180.0 




2 


30 
















2009-02-23 


180.0 




2 


20 
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Source 


Region 


ff(2000) 


5(2000) 


Date 


DIT (sec) 


K(mag) 


nod.cy 


S/N* 


SED 


reference 


alias 


SVS 20 S 


Serpens 


18 29 57.70 


+01 14 07.0 


2008-05-09 


30.0 


7.1 


2 


20 


FS 


Hai06 












2008-05-22 


30.0 




2 


20 
















2008-05-26 


30.0 




2 


25 
















2008-08-21 


30.0 




2 


30 








EC 95 


Serpens 


18 29 57.80 


+01 12 52.0 


2008-05-09 


120.0 


9.8 


2 


25 


Class II 


Hai06 












2008-05-22 


120.0 




2 


20 
















2008-05-26 


120.0 




2 


30 
















2008-08-21 


120.0 




2 


20 








Radial velocity standards 


HD 129642 




14 45 09.74 


-49 54 58.61 


2008-04-28 


120.0 


6.2 


3 


110 
















2008-05-09 


120.0 




3 


120 
















18 aug 2008 


120.0 




3 


130 
















2008-09-03 


120.0 




3 


120 








HD 105671 




12 09 54.98 


-46 12 30.20 


2008-12-28 


60.0 


5.8 


2 


150 








Gl 406 




10 56 28.86 


+07 00 52.77 


28 feb 2009 


60.0 


6.1 


2 


120 









20 



